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ABSTRACT 

We present the spectral analysis of duration- integrated broadband spectra (in ~ 30 keV— 200 MeV) 
of 15 bright BATSE gamma-ray bursts (GRBs). Some GRB spectra are very hard, with their spectral 
peak energies being above the BATSE LAD passband limit of —2 MeV. In such cases, their high- 
energy spectral parameters (peak energy and high-energy power-law indices) cannot be adequately 
constrained by BATSE LAD data alone. A few dozen bright BATSE GRBs were also observed with 
EGRET's calorimeter, TASC, in multi-MeV energy band, with a large effective area and fine energy 
resolution. Combining the BATSE and TASC data, therefore, affords spectra that span four decades of 
energy (30 keV— 200 MeV), allowing for a broadband spectral analysis with good statistics. Studying 
such broadband high-energy spectra of GRB prompt emission is crucial, as they provide key clues to 
understanding its gamma-ray emission mechanism. Among the 15 GRB spectra, we found two cases 
with a significant high-energy excess, and another case with a extremely high peak energy (E pea ^ > 
170 MeV). There have been very limited number of GRBs observed at MeV energies and above, and 
only a few instruments have been capable of observing GRBs in this energy band with such high 
sensitivity. Thus, our analysis results presented here should also help predict GRB observations with 
current and future high-energy instruments such as AGILE and GLAST, as well as with ground-based 
very-high-energy telescopes. 

Subject headings: gamma rays: bursts - gamma rays: observations 



1. INTRODUCTION 

High-energy observations of gamma-ray burst (GRB) 
prompt emission with various detectors have indicated 
that GRB co ntinuum spec t ra extend up to MeV- 
GeV energies (IHanlon et all fl994t iHurlev et all 11994 
Catelli et al.l Il998t iDingtis et all 119981 : lHarris fc Share] 
19981 iKippen et all 119981 ) GRB spectra in the MeV 
to GeV range are usually well-described by a single 
power-law with an index in the approximate range 
of -1 to -4 (ISchneid et al.l Il992t iKwok et all 119931: 
Hurlev et all Il994t IHanlon et all Il994t ISchneid et all 
1995 ][Ca"telli et al.lll998tlKippen et alJl998tlBriggs et all 
1999; IWren et ah This range is in agreement with 

the distributions of the high-energy power-law indices 
observed with BA TSE Large Area Det ectors (LADs) in 
-30 keV-2 MeV (jKaneko et al.|[200l K06 hereafter). 
However, due to the power-law nature of GRB spectra, 
photon counts above ~1 MeV are usually very low, and 
this, combined with the fact that the field of view of a 
high-energy detector is generally limited, results in much 
fewer GRBs observed in multi-MeV band than keV-band 
observations. 

The Burst and Transient Source Experiment (BATSE) 
on board the Compton Gamma-Ray Observatory 
( CGRO) has provided the largest GRB spectral database 
to date in the passband of sub-MeV range. Previously, 
K06 showed that there are some BATSE GRB spectra 
whose peak energy (Speak) lies very close to or above the 
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upper energy bound (~2 MeV) of BATSE LADs. For 
such cases, the LAD data alone cannot adequately de- 
termine either -E p0 ak or the high-energy power law index. 
Moreover, the ability to identify the high-energy power 
law component using LAD data alone is limited, since the 
LAD sensitivity d ecreases significantly toward the up- 
per energy bound (Fishma n et al.|[i989f ). Therefore, ob- 
servations of spectra extending to much higher energies 
with reasonable sensitivity are needed for these spectral 
parameters to be well determined. Combining BATSE 
LAD data with multi-MeV observations by another high- 
energy detector enables such a broadband study. 

Broadband spectral analyses of a few BATSE GRBs 
have been presented in the literatur e, using the data ob- 
tained with the CG RO instruments (|Schaefer et al.|[l998t 
iBriggs et al.lll999T ). However, those broadband spectra 
were superpositions of the deconvolved photon spectra 
that were obtained by analyzing each dataset separately 
with various photon models. Deconvolved photon counts 
are model dependent, and a spectrum constructed by 
combining individually deconvolved spectra can be quite 
different from that obtained properly by simultaneously 
fitting a common model to all datasets. 

The Total Absorption Shower Counter (TASC), the 
calorimeter of another experiment also aboard the 
CGRO - the Energetic Gamma-Ray Experiment Tele- 
scope (EGRET) spark chamber - was one of the few 
instruments capable of observing GRBs in the multi- 
MeV energy band w ith excellent energy resolution 
([Thompson et al.1 fl993). Due to its sensitivity to all di- 
rections, some bright, BATSE-triggered GRBs were also 
observed with the TASC apart from the EGRET spark 
chamber events. Since the TASC provided spectra in the 
range —1—200 MeV, broadband GRB spectra spanning 
four decades of energy can be obtained by combining 
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LAD and TASC data. Such spectra, in at least one case, 
resulted in the discovery of a distinct multi-MeV spec- 
tral component ap art from the extrapo lated sub-MeV 
BATSE component dGonzalez et al.l l2003). In this paper, 
we present duration-integrated broadband spectral anal- 
ysis of 15 bright BATSE GRBs, using LAD and TASC 
data. We first describe the instuments and available 
data types of TASC in |J2J and the selection and anal- 
ysis methodology of our study in SJ3] & |2J Then in [JSJ 
we present the results, and finally discuss the results in 

m 

2. TOTAL ABSORPTION SHOWER COUNTER (TASC) 

EGRET spark chamber was designed to observe high- 
energy gamma rays much above the BATSE energy band, 
between ~20 MeV— 30 GeV. It was equipped with an 
anti-coincidence counter and a calorimeter, TASC, which 
was located at the bottom of the module. Although the 
field of view of the EGRET spark chamber was limited 
(|Dingus et al.|[T998h . the TASC was capable of accumu- 
lating data for BATSE-triggered GRBs from all direc- 
tions, independently from the spark chamber events. 

Like the BATSE detectors, the TASC was also a 
Nal(Tl) scintillation detector but with much larger di- 
mensions of 76 cm x 76 cm with a 20-cm thickness (cor- 
responding to 8 radiation lengths). In its low-energy 
mode in the energy range of ~1— 200 MeV, the TASC 
continuously accumulated the non-burst Solar spectra 
(SOLAR), every 32.768 s. In addition, it collected Burst 
spectra (BURST) initiated by BATSE triggers, in four 
commandable time intervals of 1, 2, 4, and 16 (or 32) s. 
Both the SOLAR and BURST data provided spectra 
with 256 energy channels. An example of TASC raw 
count-rate spectra (i.e., including background counts) of 
GRB 910503 integrated over the burst duration is shown 
in Figure [T] In the spectra, the 40 K line at 1.46 MeV 
and the Fe neutron capture line at 7.64 MeV are always 
present, and they were used for on-board calibration pur- 
poses. The instrumental artifac t at ~12 MeV is due t o 
an error in the electronics design (|Thompson et all l993). 
There is also a bump around 100—200 MeV caused by 
cosmic-ray protons that passed through the TASC along 
the z-axis and deposited energy of ~ 105 MeV. The pro- 
ton spectral feature was also used to monitor the gain of 
the TASC. The FWHM energy resolution of the TASC 
is about 20% over its entire energy range. 

The response is highly dependent on the incident di- 
rection of the event photons, because of the block shape 
of the TASC Nal crystal, as well as the presence of in- 
tervening spacecraft materials surrounding the detector. 
The TASC was not capable of localizing events; there- 
fore, for GRB observations, the locations determined by 
BATSE were used to obtain detector response for each 
event. The response was calc ulated using EGS4 Monte 
Carlo code (|Nelson et alj|1985l ) with the complete CGRO 
mass model. It should be noted that the deadtime of the 
TASC is extremely high, ^60% on average. 

3. SELECTION METHODOLOGY 

We searched the TASC data for the trigger times 
of 43 bright BATSE GRBs. The 43 GRBs were se- 
lected based on the criteria of peak photon flux > 
10 photons cm~ 2 s _1 in the 1024 ms timescale and 
BATSE channel-4 energy fluence (> 300 keV) > 5.0 x 



10~ 5 ergs cm -2 . The flux and fluenc e values were 
taken from the 4B (|Paciesas et al.1 1 1999f ) and the cur- 
rent 5 BATSE catalog. Among these, we identified 28 
GRBs to have significant detections in the TASC data. 
In this analysis our sample consists of 15 bursts (out of 
28), for which we had computed TASC detector response 
matrices prior to this study. The 15 GRBs are listed in 
Table [1] along with the time and energy intervals used for 
analysis, for which the selection methods are described 
below. 

For this analysis we used TASC SOLAR data, which 
provides 32.768-s time resolution. Since all of the 15 
events are very strong and most of them have durations 
longer than 33 seconds, the time resolution is adequate 
for the duration-integrated spectral analysis performed 
here. As for the selection of corresponding BATSE data, 
we use the LAD Continuous (CONT) data of the bright- 
est LAD, with 2.048-s time resolution. The brightest 
LAD (i.e., the LAD that recorded the highest counts) 
for each event is also listed in Table □ The CONT 
data were able to provide time intervals with a sufficient 
match to the SOLAR data time intervals, which often 
began before the BATSE trigger time. The time inter- 
vals chosen for analysis were determined based upon the 
detection significance above background for either the 
LAD or TASC lightcurves. For all but one event (with 
a very weak tail in the LAD data), the time intervals 
chosen for the analysis contained their BATSE Tgo in- 
tervals. All 15 events are also included in the BATSE 
LAD spectral analysis presented in K06, although the 
data type used here (hence the time and the energy in- 
tervals) may be different. In Figure[2l we show, for all 15 
GRBs, both the LAD and TASC lightcurves over their 
entire energy ranges. In the Figure, the background mod- 
els and the time intervals used for the analysis are also 
shown. The background models were determined by fit- 
ting a low-order polynomial function to the spectra accu- 
mulated for several hundreds of seconds before and after 
the burst duration interval. 

LAD CONT data provided 16 energy channels in the 
energy range of ~30 keV— 2 MeV. The lowest few chan- 
nels are usually below the electronic threshold, and the 
highest channel is an energy overflow channel; therefore, 
a total of about 13 energy channels were usually included 
in the analysis. For the TASC SOLAR data, the lowest 
6—7 channels are always excluded to assure the exclusion 
of an electronic cutoff. This translates into a lower energy 
bound of ^1.3 MeV (TableHJ column 9). In addition, the 
uppermost 10—20 channels are also excluded, depending 
upon the gain of the detector at the time of the event. 
The resulting upper bound for the TASC energy range 
was ~130-200 MeV (Tabled! column 10). Notice that 
overlap in energy between the LAD and TASC datasets 
of a few hundred keV exists in each event. 

4. SPECTRAL ANALYSIS 

We converted the TASC data into the BATSE BFITS 
format, and used the BATSE spectral analysis tool 
RMFIT 6 for this analysis. For comparison purpose 
and consistency, the same set of photon models used 

5 See |http77/ gammaray,nsstc. nasa.gov/batse/grb/catalog/c urrenty] 

6 R.S. Mallozzi, R.D. Preece, & M.S. Briggs, " RMFIT, ~K 
Lightcurve and Spectral Analysis Tool," © Robert D. Preece, Uni- 
versity of Alabama in Huntsville 
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to analyze the bright BATSE GRBs in K06, namely, 
a single power-law (PWRL), empirical GRB functions 
with and without the hig h-energy power-la w (BAND 
and COMP, respectively; iBand et all 11993ft. and the 
smoothly-broken power-law ( SBPL; IPreece etallUQfll 
lRvddll999t IPreece et ai]l2000| ) with its subsets, were fit- 
ted to the joint LAD-TASC dataset. Extensive discus- 
sions on these spectral models are found in K06. Since 
we are only analyzing the duration-integrated spectra, 
the GRB function with fixed high-energy index (3 (BETA 
model in K06) was not applicable and thus not used. 
The PWRL and COMP models are expected to result in 
very poor fits to the joint spectra because of the broad 
energy coverage, as well as the fact that our sample con- 
tains the very brightest of all BATSE GRBs, which have 
been shown to have a smoothly-broken power-law behav- 
ior with no cutoff in keV band (their best-fit models in 
the LAD-only analysis are SBPL or BAND; K06). 

In order to account for uncertainties in the effective 
areas of the two detectors, we added an effective area 
correction (EAC) term to the spectral model of each fit: 
this is a multiplicative factor, to normalize the data of 
the second detector to those of the first detector. In our 
case, the EAC factor always normalized the TASC data 
to the LAD data. The EAC values can vary from event 
to event and were determined by simultaneously fitting 
the BAND model to both datasets with the EAC term 
as a free parameter. Once the EAC factor was found, it 
was kept fixed in the subsequent final fits. Between the 
LAD and the TASC datasets, they were always found 
to be unusually large (~0.5). We investigate this issue 
further in £15.11 

5. ANALYSIS RESULTS 

As we predicted, broadband high-statistics data confi- 
dently rejected PWRL and COMP as best-fit models for 
all 15 spectra in our sample. Following the LAD GRB 
analysis in K 06, wc determined the best-fit (BEST) mod- 
els for the 15 joint spectra as well. The BEST model is 
the simplest, statistically well- fit model among the four 
spectral models described above that is required by the 
observed spectrum. Since we have four models with 2, 3, 
4, and 5 free parameters, we compared resulting \ 2 val- 
ues between all combination of models, and searched for 
significant (>99.9%) improvement in y 2 as we move from 
the simplest to more complicated models. We also take 
into account the spectral parameter constraints. The de- 
tails of determining the BEST models is found in K06. 

For all GRBs, the joint count spectra and the BEST 
models, along with the corresponding vF v spectra are 
shown in Figure [3l The spectral parameters of the 
BEST model arc also presented in Table [2l In the ta- 
ble, £ p k value 7 is the peak energy in vF v spectrum and 
break energy (Eh) value is the break energy of a broken 
power law, regardless of the model fitted. As mentioned 
in K06, -Epcak and E\, are not necessarily the same for 
a given spectrum due to the curvature around the break 
energy. The derivation of these values are presented in 
appendices of K 06. It must be noted that due to the high 
photon counts of these events, uncertainties in the data 

7 Hereafter, we use £ pk notation to represent the actual peak en- 
ergy in vF v spectrum, to distinguish from a fit parameter -E pe ak m 
BAND model. E paak = £ pk only if p < -2. 



are dominated by systematics, which can be large espe- 
cially at lower energies (< 100 keV) in the LAD data. 
This can result in a relatively large \ 2 value even for an 
acceptable fit, which can be seen in the sigma residuals 
in Figure [3] 

To compare these parameters with the parameter dis- 
tributions of the larger sample of bright BATSE bursts 
presented in K06, we overplot in Figures [4] and [5] the 
BEST model parameters of the jointly-analyzed events 
on top of the time-integrated LAD spectral parameter 
distributions, taken from K06. It is evident in terms of 
the photon fluence (in 25—2500 keV; FigurelH top panel) 
that these 15 events are in the very brightest group of 
all BATSE GRBs. No bias or tendency is seen in the 
distributions of spectral indices (Figure [4l bottom pan- 
els). It is clear, however, that the £ p k for the 15 events 
belong to the higher end of the BATSE distributions 
(Figure [5]), and by themselves forms a quasi-Gaussian 
with a median of (£ p k) = 517 keV. A likely reason for 
this is because our sample selection was based on the 
high photon flux and fluence above 300 keV: a group 
correlation between burst brightness and £ p k was previ- 
ously identified using an early BATSE GRB sample by 
Mallo zzi et al.l (|1995h . We note also that because we se- 
lected only GRBs with significant detections in TASC 
data (> 1 MeV), this naturally introduces a preference 
for GRBs with higher £ p k to be included in our sample. 
We found that the photon fluence and the energy fluence 
determined with the joint spectra were consistent with 
the values estimated by extrapolating the LAD-only fit 
spectra up to 200 MeV. In addition, we did not find sig- 
nificant correlation between photon fluence and .Epoak or 
E\, within our sample. 

To illustrate the improvements in parameter con- 
straints as a result of the joint analysis, in comparison 
with the single-detector analysis, we show the spectral 
parameters determined by the joint fits and the individ- 
ual detector fits in Figure [(51 In the single-detector anal- 
ysis, the LAD data were fitted with the BEST model 
listed in Table El while the TASC data were fitted with 
PWRL with pivot energy of 10 MeV. The PWRL indices 
of TASC are compared with the high-energy indices of 
the BEST models in Figure [6] (top right panel) . It is not 
surprising that in most cases the spectral parameters are 
determined much better with the joint analysis than with 
the individual cases. We also note that spectral parame- 
ters of our TASC-only fits were consistent with the pre- 
vious TASC analysis results found in the literature (e.g. , 
Schneid et alJll992HKwok et al.|1993tlHurlev et al|1994 



Catelli et afl 119981 : iBriggs et all [19991 : IWren et all 12002 
within a few sigma uncertainties. Moreover, a few of 
the 15 GRBs were also detected by the EGRET spark 
chamber in even higher energies, w ith their flux values 
consistent with the TASC spectra ( Schneid et al.l 119921 : 
iKwok et al.lll993l : iHurlev et al.l[l99l ~ 

As seen in Figure [6] (top left panel) , the low-energy 
spectral indices are already well constrained by the LAD 
data alone and are least affected by the addition of TASC 
data. On the other hand, the high-energy indices are 
much better constrained by the joint analysis, as one 
would expect. The values determined by the joint anal- 
ysis are nearly always found between the values derived 
from the LAD-only and TASC-only analyses. In cases 
where the LAD-determined index differs from the jointly- 
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determined index by more than a few a (i.e., triggers 249 
and 2831), the spectral break energy was also found to 
change significantly (Figure [51 bottom right panel). The 
£ p k values found by the joint analysis are consistent with 
those found by the LAD-only fits, although in many cases 
it seems to settle in the higher ends of the value ranges 
determined by LAD. We observed one case, trigger num- 
ber 3523, in which the value of £ p k determined with the 
joint fit was less constrained than the one found with the 
LAD-only fit. This is due to the high-energy index being 
very close to —2 (—2.01 ± 0.03), making the constraint 
of £ p k very difficult, by definition (see its vF v spectrum 
in Figure [3|). There was one event (trigger number 2329) 
for which the £ p k could not be determined even with the 
joint spectra because the high-energy index was above —2 
by 7.5er (see Table and therefore, the vF v spectrum 
did not peak within the energy range. This may indicate 
a lower limit in £ p k of 167 MeV (upper energy bound 
of the spectrum) for this event. The lightcurve of this 
burst (Figure [2|) shows that in the first time interval (T- 
33 to s), the TASC flux is clearly dominant and thus, 
spectrally harder than the second interval spectrum. We 
note that the time-resolved LAD-TASC joint analysis of 
this event in a time interval of 1—23 seconds also found 
that the high-energy indices were always above —2 by at 
least la (IGonzalez et al.ll2004t iKaneko et aLll2005[ K 

Although the spectral fits presented in Table [2] were 
all sufficiently good (large x 2 were due to systemat- 
ics in LAD), indications of a high-energy excess were 
found in the residual patterns of at least three spectra 
(GRBs 920902, 941017, and 980923; triggers 1886, 3245, 
and 7113). This was much more evident in the spec- 
trum of 3245: The high-energy excess above 10 MeV 
is seen in sigma residuals of the spectrum (Figure [3]), 
and this is also evident from its lightcurves (Figure . 
For this event, fitting an additional high-energy PWRL 
together with the BEST model resulted in an improve- 
ment in x 2 °f 18-1 (for Adof = 2), corresponding to a 
chance probability of 10~ 4 . In the case of 7113, the x 2 
improvement for adding high-energy PWRL was 14.0, in- 
dicating a chance probability of ~10 -3 , whereas for the 
last event (1886), the corresponding x 2 improvement was 
only about 7, which translates into a chance probability 
of 0.03. Although, in all case, the additional high-energy 
PWRL indices could not be well determined (la uncer- 
tainties > 0.5), the indices were rather hard, likely > — 2. 
Time-resolved analysis of 3245 previously revealed, with 
much higher significance, a high-energy spectral compo- 
nent th at deviates from the e xtrapolated keV LAD com- 
ponent (jGonzalez et alJl2003f >. This high -energy compo- 
nent emerged later than and remained bright longer than 
the keV spectral component. Although our analysis pre- 
sented here is only for duration-integrated spectra, the 
delayed nature of the excess high-energy emission is clear 
from the lightcurves of all three events (Figure [2|) . 

5.1. Effective Area Correction Issue 

All detector response models have uncertainties. To 
allow for the uncertainties in their effective areas, we 
employ normalization factors between detectors when si- 
multaneously analyzing spectra from multiple detectors. 
Usually o nly about 10% dif ference between datasets is 
expected (jBriggs et al.ll 1999ft . In our joint analysis, how- 
ever, we find the discrepancy between the LAD and 



TASC data to be relatively large, ^30—80%, as indi- 
cated by the EAC factors (see Table EJ). Disagreements 
between LAD and TASC data were also previously found 
in some of the co mposed spectra (jSchaefer et al.l [l998t 
iBriggs et al.lll999i ). although the response matrices used 
here have been newly calculated. 

As mentioned earlier, the TASC effective area de- 
pends highly upon the incident angles of events, be- 
cause of the geometric area of the Nal crystal as well 
as varying amount of intervening spacecraft material 
(jThompson et al.1 [T993). Therefore, we investigated the 
EAC factors in our sample in terms of the incident an- 
gles, as well as other potential contributing factors, such 
as the TASC live time, event brightness, energy range, 
and spectral parameters. We found no apparent correla- 
tion in any of those; however, we noticed a striking re- 
semblance between the plots of EAC vs. incident zenith 
angle and of the effective area vs. zenith angle. The 
comparison is shown in Figure in which our EAC val- 
ues for the 15 GRBs and the normalized TASC effective 
area are plotted against incident zenith angles. This in- 
dicates that the discrepancy between the two detectors 
is more severe at the zenith angle where the TASC ef- 
fective area is smaller. As a matter of fact, a disagree- 
ment between the calculated effective area and the ac- 
tual experimental value was found at the time of the 
TASC instrument calibration, which was attributed to 
the CGR O mass model underest imating the intervening 
material (jThompson et al.l fl993h . The EAC factors in 
our analysis were always found to be less than 1, mean- 
ing count rates in the TASC data are overestimated. The 
count overestimation becomes more apparent when there 
is larger amount of intervening material, namely, when 
the effective area is smaller. Consequently, we conclude 
that the CGRO mass model used to determine the effec- 
tive area indeed underestimates the intervening material, 
at each zenith angle, resulting in overestimation of the 
photon counts observed with TASC. 

6. SUMMARY AND DISCUSSION 

In this study, we extended the BATSE GRB spec- 
tral analysis (K06) to high-energy broadband spectra 
obtained by combining BATSE LAD data and EGRET 
TASC data. Time-integrated joint spectra of 15 hard 
BATSE GRBs were analyzed in order to probe high- 
energy spectral properties of prompt emission. The 
TASC data in multi-MeV energy band with fine en- 
ergy resolution, clearly confirmed that the GRB spec- 
tra do extend up to ^200 MeV and probably beyond. 
The joint broadband spectra in the energy range of 
~30 keV— 200 MeV indeed constrain the high-energy 
spectral indices and break energies of strong GRBs that 
have significant MeV emission much better than single 
detector analysis. In most cases, £ p k (and E\f) values 
derived with LAD data alone and values found by joint 
analysis were consistent within la uncertainty. However, 
in a few cases, the jointly-fit indices differed significantly 
from those determined with the LAD data alone. This 
indicates the possibility that some high-energy indices 
obtained with BATSE LADs alone may not reflect the 
intrinsic values, and thus highlights the importance of 
broadband spectral analysis. 

We identified one case (GRB 930506; trigger 2329) in 
which the £ p k could be extremely high, with a lower limit 
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of ~167 MeV. We note here, however, that this possibly 
very high £ p k should be interpreted with caution for two 
reasons: 

1. While the high-energy power law index being > — 2 
is statistically very significant (by 7.5c) as deter- 
mined by the statistical error bar, there is the possi- 
bility of a systematic error existing in a joint fit be- 
tween two instruments, which is not taken into ac- 
count here. Including the systematic uncertainties 
in the analysis could bring down the index value 
to < —2, in which case the £ p k would probably 
become a few MeV. 

2. Another possibility is that there are two spectral 
components in this spectrum; the "regular" one 
with £ p k of a few MeV and an additional high- 
energy multi-MeV component. In addition, if the 
spectra evolved similar to the case of GRB 941017, 
the duration-integrated spectrum could be a su- 
perposition of two evolving components. In such 
cases, the single model used here could hinder cor- 
rect identification of the £ p k as well as the high- 
energy power law index. We found, however, no 
statistical evidence for such high-energy component 
in the spectrum of GRB 930506. 

In any case, the £ p k of GRB 930506 is clearly above the 
BATSE LAD energy range. This, combined with the fact 
that in some cases the £ p k values were found to be slightly 
higher than those determined by LAD data only, may 
indicate that there exists a tail population of £ p k values 
extending to a few hundred MeV, especially given the 
limited sensitivity of TASC. 

In addition, while the broadband spectra were mostly 
consistent with broken power laws (BAND or SBPL), 
indications of high-energy excess were also found in two 
events (GRBs 941017 & 980923; triggers 3245 & 7113). 
For one of them (GRB 941017; trigger 3245), the dis- 
tinct MeV component has been clearly identified wit h 
time- resolved spectral analysis (|Gonzalez et al.l [2003), 
with much stronger significance than with the time- 
integrated spectrum. For this particular GRB, the ex- 
istence of the extra component was further reinforced by 
the COMPTEL spe ctra, which covered an energy range 
of 30keV-10MeV (jKaneko et al.ll200l . 

Possible explanations that have been proposed for 
such a high-energy component include the Comp- 
ton upscattering of synchrotron photons in the re- 
verse shoc k by the synchrotron- emitting relativistic 
electrons (|Granot fc Guettal l2003h . and the electro- 
magnetic cascade emission of ultra-relativistic baryons 
through photo-pion interactions and subsequent pion de- 

8 http:/ /agile. iasf-roma.inaf.it/ 

9 http://wwwmagic.mppmu.mpg.de/ 

10 http:/ /veritas. sao.arizona.edu/ 

11 http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html 

12 http://glast.gsfc.nasa.gov/ 



cay dDermer fc Atovanll2004f) . In the synchrotron shock 
model |Katdll994l ). the synchrotron self-Compton emis- 
sion associated with the synchrotron component that 
peaks in the keV band is expected to be observed in 
the MeV to TeV energy band, depending on the Lorentz 
factor of the electrons (jGuetta fc Granotll2003l ). How- 
ever, the brightness and delayed nature of the high- 
energy component are inconsistent with the synchrotron- 
self Compton. It is particularly interesting if the compo- 
nent is indeed due to the relativistic baryons, since the 
observed component may be direct evidence of baryonic 
acceleration, namely, cosmic rays. 

The high-ener gy power-law compone nt observed in 
GRB 941017 by fconzalez etahl (|2003D was described 
with an additional power law of index ~ —1. This re- 
quires that there exists another break energy (and a true 
£ p k) above 200 MeV, in order to avoid the energy diver- 
gence. Although the redshift for this event is unknown, 
a rough upper limit energy for such a break could be 
placed at ~1 TeV, solely by the total (isotropic equiva- 
lent) energy constraint of ~ 10 53 ergs, and by assuming 
the burst originated nearby (z <C 1). 

As was the case for GRB 941017, time-resolved spectral 
analyses of high-energy GRBs are needed in order to ex- 
plicitly identify similar distinct spectral component. We 
have indeed performed such time-resolved analyses for all 
GRBs presented here, the result of which is the subject of 
another paper that will follow this work (M.M. Gonzalez 
et al. in preparation). It is possible that the time- 
resolved spectral analysis of broadband spectra would re- 
veal such high-energy spectral components in many more 
GRBs. Also in the MeV energy band and above, there 
should be a spectral cu toff due to pair-production at- 
tenuation (Baring 2000). Such a cutoff energy would 
determine the bulk Lorentz factor T, which is one of the 
main factors, along with a redshift, preventing the deter- 
mination of the int rinsic £ p k values in the r est frame of 
expanding matter (Zhang Sz Meszaros 2004). 

Unfortunately, the TASC data are only available for 
the brightest of BATSE GRBs due to its sensitivity lim- 
itation. Although GRB spectra with a distinct MeV 
component seem to be rare, and we observed no pair- 
production cutoff in our analysis, these high-energy spec- 
tral features could be sought after using currently ex- 
isting very high-energy telescopes, such as AGILE 8 , 
MAGIC 9 , VERITAS 10 , and HESS 11 . Upcoming obser- 
vations by GLAST 12 with much higher sensitivity in an 
unprecedented broad energy band of 10 keV— ~30 GeV, 
are anticipated to reveal broadband spectral character- 
istics of GRBs, which holds significant clues to under- 
standing GRB emission mechanism. 
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Fig. 1. — An example TASC raw count spectrum of GRB 910503, for the time interval including the burst episode (—9.8 to 121 s since 
BATSE trigger time). The count rate includes background and the source counts. 40 K line, Fe line, the instrumental artifact, and the 
cosmic-ray proton features are clearly seen. 



TABLE 1 

Time and energy intervals used for. the LAD-TASC joint spectral analysis of 15 GRBs. Time intervals are shown in units 

OF SECOND SINCE BATSE TRIGGER. 









Time Interval 


LAD Energy 


TASC Energy 


GRB 


BATSE 


LAD 


Start 


End 


Start 


End 


Start 


End 


Name 


Trigger No. 


No. 


(s) 


(s) 


(keV) (MeV) 


(MeV) (MeV) 


910503 


143 


6 


-9.8 


88.5 


36.6 


1.81 


1.29 


201 


910601 


249 


2 


15.9 


48.6 


37.0 


1.82 


1.32 


198 


920525 


1625 


5 


-18.4 


47.1 


32.6 


1.85 


1.30 


195 


920622 


1663 


4 


-30.9 


34.6 


31.9 


1.81 


1.32 


195 


920902 


1886 


5 


-10.2 


120.8 


32.7 


1.85 


1.39 


196 


930506 


2329 


3 


-33.1 


32.5 


44.2 


1.85 


1.18 


167 


940217 


2831 


2 


-8.2 


188.4 


37.1 


1.82 


1.25 


187 


940703 


3057 


1 


-1.0 


97.2 


33.3 


1.88 


1.28 


139 


940921 


3178 


2 


-25.6 


72.6 


28.1 


1.82 


1.31 


178 


941017 


3245 


4 


-51.2 


210.9 


31.8 


1.81 


1.25 


199 


950425 


3523 


6 


-23.6 


107.5 


36.3 


1.80 


1.57 


197 


970315 


6124 


2 


-15.8 


82.5 


37.1 


1.82 


1.51 


189 


980923 


7113 


7 


-20.0 


78.7 


33.4 


1.92 


1.44 


131 


990104 


7301 


7 


-4.5 


257.7 


29.3 


1.96 


1.45 


202 


990123 


7343 





-0.1 


98.3 


33.1 


1.81 


1.40 


128 
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Fig. 2. — LAD and TASC lightcurves of all 15 GRBs. The dotted lines show background models and the vertical dashed lines indicate 
the selected integration time intervals. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 3. — Count spectra and the BEST models with sigma residuals (left panels) and corresponding vF v spectra (right panels) of LAD- 
TASC joint fits of all 15 GRBs. The BEST models are shown in solid lines in both panels. The TASC data have been binned for display 
purpose only. 
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Fig. 3. — Continued. 
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Fig. 3. — Continued. 
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Fig. 3. — Continued. 
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Fig. 4. — Parameter distribution comparisons of 15 joint events with 342 bright BATSE GRBs. The crosses represent photon fluence 
(top) and spectral indices (bottom) of jointly analyzed events, and histogram shows the distribution of LAD-only analysis events. 
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Fig. 5. — Parameter distribution comparisons of 15 joint events with 342 bright BATSE GRBs. The crosses represent £p k (left) and 
E h (right) values of jointly analyzed events, and histogram shows the distribution of LAD events. In the £ pk distribution, GRB 930506 
(trigger 2329) is excluded due to a very high £ pk > 167 MeV. 
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TABLE 2 

Summary of LAD-TASC joint fit results for time-integrated spectra. Full functional form and parameter description of 

EACH SPECTRAL MODEL CAN BE FOUND IN K06. lcr UNCERTAINTIES ARE SHOWN IN PARENTHESES. 



Spectral Fit Parameters 



GRB 


BATSE 


BEST 




Amplitude 


1 


2 

Low Index 


Hi 


;5 

?h Index 


4 


Break 


EAC 


X 2 /dof 


Name 


Trig No. 


Model 


(ph s 


1 cm" 2 ) 


(kcV) 














(keV) 


Scale, A 








910503 


143 


SBPL 





0116 


(0.0000) 


612 ( 28) 


-1 


18 


(0.01) 


-2 


25 


(0.03) 


467 ( 17) 


0.20 





63 


204.4/224 


910601 


240 


BAND 





0567 


(0.0003) 


515 ( 8) 


-1 


04 


(0.01) 


-3 


08 


(0.10) 


562 ( 26) 







71 


288.0/224 


920525 


1625 


BAND 





0151 


(0.0003) 


427 ( 18) 


-0 


88 


(0.03) 


-2 


45 


(0.10) 


311 ( 23) 







19 


250.1/221 


920622 


1663 


SBPL 





0174 


(0.0001) 


386 ( 20) 


-1 


11 


(0.01) 


-2 


24 


(0.04) 


286 ( 11) 


0.20 





73 


228.9/220 


920902 


1886 


BAND 





0043 


(0.0001) 


469 ( 26) 


-0 


41 


(0.05) 


-2 


25 


(0.11) 


284 ( 21) 







70 


110.8/219 


930506 


2329 


SBPL 





0220 


(0.0001) 


> 167000? 


-1 


26 


(0.01) 


-1 


85 


(0.02) 


349 ( 22) 


0.30 





50 


226.6/208 


940217 


2831 


BAND 





0203 


(0.0001) 


553 ( 12) 


-1 


17 


(0.01) 


-2 


88 


(0.10) 


580 ( 35) 







49 


257.7/222 


940703 


3057 


SBPL 





0315 


(0.0001) 


585 ( 29) 


-1 


15 


(0.01) 


-2 


27 


(0.03) 


393 ( 13) 


0.30 





33 


211.8/206 


940921 


3178 


BAND 





0067 


(0.0002) 


684 ( 65) 


-1 


17 


(0.03) 


-2 


30 


(0.13) 


478 ( 74) 







73 


218.6/217 


941017 


3245 


SBPL 





0086 


(0.0001) 


250 ( 17) 


-1 


07 


(0.02) 


-2 


07 


(0.03) 


185 ( 7) 


0.10 





44 


183.4/215 


950425 


3523 


SBPL 





0095 


(0.0001) 


2150 (1568) 


-1 


03 


(0.02) 


-2 


04 


(0.03) 


479 ( 30) 


0.30 





69 


203.8/215 


970315 


6124 


BAND 





0143 


(0.0004) 


392 ( 19) 


-0 


90 


(0.03) 


-2 


77 


(0.29) 


345 ( 55) 







63 


207.3/215 


980923 


7113 


BAND 





0639 


(0.0005) 


379 ( 5) 


-0 


66 


(0.01) 


-2 


65 


(0.05) 


294 ( 8) 







37 


267.7/198 


990104 


7301 


BAND 





0190 


(0.0006) 


696 (119) 


-1 


31 


(0.05) 


-2 


40 


(0.09) 


560 (104) 







35 


202.5/217 


990123 


7343 


SBPL 





0256 


(0.0001) 


517 ( 16) 


1 


04 


(0.01) 


-2 


58 


(0.03) 


436 ( 10) 


0.30 





68 


400.4/200 



1 Energy at which the vF v spectrum peaks. For BAND, it is a Fitted, E pc!i k, if /3 < —2, and for SBPL, it is a calculated 
£ P cak (see Appendix B in K06). 

2 a for BAND and Ai for SBPL. 

3 13 for BAND and A 2 for SBPL. 

4 Fitted E h for SBPL, and calculated E h for BAND (see Appendix C in K06). 



